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Basing on self-consistent solution of non-linear GL-equations, the phase boundary is found, which 
divides the regions of I- and II-order phase transitions of a superconducting cylinder in magnetic 
field to normal state. This boundary is a complicated function of the parameters (m, R, x) (m is 
the vorticity, R is the cylinder radius, x is the GL-parameter) , which does not coincide with the 
simple phase boundary x = 1/ v2, dividing the regions of I- and II-order phase transitions in infinite 
(open) superconducting systems. 



I. INTRODUCTION 

The GL-theory is widely used for studying the gen- 
eral properties of the superconducting state. This the- 
ory leads to two coupled non-linear equations for the or- 
der parameter tp and the magnetic field vector-potential 
A, which are usually solved using various simplifying as- 
sumptions. In a number of papers [2-8] the particular 
case was considered of a long superconducting cylinder 
of radius R, placed in the axial magnetic field H . In 
this case the 3-dimensional GL-equations reduce to one- 
dimensional form, which enables one to find numerically 
the exact self-consistent solutions. In such way it is possi- 
ble to study specific non-linear effects, as well as the role 
of the sample boundary. For instance, it was shown in 
[6,7], that one-dimensional solution for the order parame- 
ter ip (with fixed vorticity m and varying H) may change 
its form either gradually (in one interval of parameters 
(R, x) , x is the GL-parameter) , or abruptly (in the other 
interval of (R, x)), undergoing I-order jump transforma- 
tion. Such jump transformations, in principle, may be 
observable, because they are accompanied by jumps of 
the magnetisation, M(H). 

In the present paper the phase boundary is found, 
which divides the region of parameters (R, x), where the 
superconducting solution (of fixed m) ends up (in the in- 
creasing external field) by I-order jump to normal state 
(-0 = 0), from the region (R, x) , where the solution van- 
ishes gradually, by II-order phase transition. This phase 
boundary is a complicated function of R and x, different 
from the simple boundary x = 1/ a/2, which divides the 
I- and II-order phase transitions in infinite (open) super- 
conductors [9] . Other topics are also touched on (such as 
metastability, the paramagnetic Meissner effect, the pin- 
ning of vortices to the sample boundary, the linearized 
equation approximation, etc.). 

In Sec. II the problem is formulated and the basic 
GL-equations, used in calculations, are written. Sec. II 
contains the numerical results, alongside with necessary 



comments. In Sec. IV the results are summarized and 
discussed. 



II. EQUATIONS 

Below the case is considered of a long superconduct- 
ing cylinder of radius R, in the external magnetic field 
H > 0, which is parallel to the cylinder element. In the 
cylindrical co-ordinates the system of GL-equations may 
be written in dimensionless form [6] 
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Here U(p) is the dimensionless field potential; b(p) is the 
dimensionless magnetic field; ip(p) is the normalized or- 
der parameter; p — r/A, A is the field penetration length; 
A = x£, where £ is the coherence length, x is the GL- 
parameter. The dimensioned potential A, field B and 
current j s are related to the corresponding dimensionless 
quantities by the formulae [6]: 
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(The field B in (3) is normalized by H x = (j) /(2n\ 2 ), 
with b = B/Hy, instead of H\ one can normalize by 
= 0o/(2< 2 ), or by H K - O /(2<A) = H e /x. The 
coefficients in (I), (2) would change accordingly.) The 
vorticity m in (3) specifies how many flux quanta are 
associated with the vortex, centered at the cylinder axis 
(the so-called giant-vortex state). 

The boundary conditions to Eq. (I) are [7]: 
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where p\ = R/X, h\ = H/H\. 

The boundary conditions to Eq. (2) are [7]: 
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The magnetic moment (or, magnetisation) of the cylin- 
der, related to the unity volume, may be written in a form 
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where -B at) is the mean field value inside the superconduc- 
tor, $i is the total magnetic flux, confined in the cylin- 
der. In the normalization (3), denoting b = B av /H\, 
h x = H/H x , Mx = M/Hx, one finds 



4nMx = b-hx, b= -^(t/i+m), (6) 
Pi 
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Accordingly, the normalized Gibbs free energy of the 
system may be written as [7] 
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Here AG = G s — G n is the difference of free energies 
in superconducting and normal states; 6(0) = B(Q)/Hx, 
B(0) is the magnetic field at the cylinder axis; H cm = 
0o/(27rv / 2A^) is thermodinamical critical field of massive 
superconductor; go is the condensation energy with ac- 
count for the order parameter gradient. The expressions 
(6), (7) may be used for calculating the corresponding 
quantities. [Instead of pi the notation Rx — R/X = pi 
will be used below. 1 



III. NUMERICAL RESULTS 

The solutions of Eqs. (l)-(5) depend on the space co- 
ordinate p and several parameters, for instance, tp(p) = 



tp(m, Rx, x, hx; p) (analogously for the potential U(p) 
and the field b(p)). Let the vorticity m be fixed (m — 
0, 1, 2,...) and consider at first the case m = (the 
vortex- free Meissner state). Consider the plane of pa- 
rameters (Rx,x) (see Fig. 1(a)). In every point of this 
plane there exists a set of solutions of Eqs. (l)-(5), which 
depend paramctrically on the external field hx- (Several 
points, laying along the line Rx = 4 are numerated as 
1-6.) One may envisage a peep-hole, pearced in every 
point (Rx, x), which allows to see the content of the cor- 
responding sub- volume. The set of solutions tp(hx; p) is 
unique for each sub- volume and may be characterised, for 
instance, by the field dependence of the maximal value 
of the order parameter tp max (hx), or by the form of the 
magnetisation curve Mx(hx) (6). The examples of such 
dependencies in different points of the plane (Rx, are 
given in Figs. 1(6, c) (only the case hx > is considered; 
some illustrations for the case hx < 0, as well as the cor- 
responding co-ordinate dependencies, may be found in 
[6,7])- 

It is clear from Fig. 1(6), that the characteristic be- 
havior tp m ax(hx) depends essentially on the value of k. 
For small x, the value tp max (hx) terminates (the curves 
1-4) by jump at some point hx = h s , where (if hx is 
increased) the I-order phase transition to normal state 
(tp(p) = 0) occures. The region, where the supercon- 
ducting solutions terminate by I-order phase transition, 
is marked in Fig. 1(a) as sj. 

For larger x (the curves 5,6) there is also a jump in 
tpmax(hx) at some point hx = h s , but with a "tail" ap- 
pearing on the curve. If the field hx increases, the super- 
conducting solutions 5,6 vanish gradually at the point 
h c , by Il-order phase transition to normal state. The 
region, where the superconducting solutions terminate 
by Il-order phase transition, is marked in Fig. 1(a) as 
sn- [The appearance of the tail on the magnetisation 
curve means the transition of the solution to the edge- 
suppressed form, see [7] for details.] 

It is evident, that for small radius cylinder (Rx < 
1.69) the superconducting solution terminates by Il-order 
phase transition, even in type-1 (i.e. small x) supercon- 
ductors [10]. The transformation of the solutions with 
diminishing radius Rx is illustrated in Fig. 2 for x = 0.7. 

Note, that if the line x = 1 in Fig. 1 (a) is followed from 
large to small Rx, the superconducting states, which lay 
along this line, display at first the Il-order phase tran- 
sition in magnetic field (for larger Rx), then the I-order 
(for intermediate Rx), and again the Il-order (for smaller 
Rx)- Only, if x > 1.05, all the solutions display Il-order 
behavior. 

Notice also, that the state to = is totally diamagnetic 
(-AttMx > 0). 

Because in every point of sn-region the order param- 
eter vanishes by Il-order phase transition (tp ma x — * 0, 
see Fig. 1(6)), the superconducting phase boundary in 
magnetic field, h c , may be found analytically, by lineariz- 
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ing the system (1), (2) (with account, that ip -C 1 and 
6 :=y h\), and passing to single linear equation for the 
order parameter [11] , whose solution may be expressed 
in terms of the Kummer functions (see also [3,4,5,12]). 
However, inside si-region (where the solution terminates 
by jump from a finite value ipmax to zero) the phase 
boundary h s (i.e. the highest field h\ still compatible 
with the superconductivity) can not be found by solving 
the linearized equation, but full system (l)-(5) is needed. 

The analogous investigation can be carried out in the 
case m = 1 (see Fig. 3), with a single vortex on the 
cylinder axis. 

In Fig. 3(a) are shown: the region si, where the su- 
perconducting state terminates (if the field is increased) 
by I-order jump to the normal state, having finite value 
ipmax at the transition point; the region sn, where the 
superconductivity vanishes by Il-order phase transition; 
the curve Si_n, which represents the boundary between 
I- and Il-order phase transitions. 

The behavior of the order parameter ip max {hx) and 
of the magnetisation M\(h\) in different points of the 
plane (R\, x) are shown in Figs. 3(6, c) (and in Fig. 4). 
For small x (the curves 1,2) the solutions terminate by 
I-order jump. When the line si_n is crossed, the tail 
appears on the curves 3,4, which widens, if R\ and x 
increase. If Rx diminishes (Fig. 4), the magnitude of the 
jump in ipmax also diminishes, and the solutions termi- 
nate (if the field is increased) by Il-order phase transition 
to normal state. 

On the curve C ns (Fig. 3(a)) the value i\) max = 0. The 
letter n denotes the normal metal region (ip = 0); here 
the superconducting state (to = 1) is impossible. [In this 
region the radius R\ is too small, and the vortex own 
field is too strong to be confined within the mesoscopic 
sample] It is evident, that when the radius R\ dimin- 
ishes (but x is fixed) the transition from s- to n-state 
always is Il-order phase transition, however the width of 
the region between the curves si_n and C ns (where Il- 
order transition exists) gets very small for small x. The 
curve C ns may be well approximated by the dependence 
R\ ~ a/x (or R^ = xR\ — a), with a = 1.34. 

Notice, that in any point of s-rcgion in Fig. 3(a) 
the magnetisation function M\(h\) (Fig. 3(c)) has two 
parts: the paramagnetic (M\ > 0) and diamagnetic 
(M\ < 0). This is because the superconducting cur- 
rent has two components, j s = j p + jj. One of these 
currents (j p ) screens the own field of the vortex (to = 1) 
and flows around the vortex axis in counter-clockwise di- 
rection (the paramagnetic current). The second current 
(jd) screens out the external field h\ and flows near the 
cylinder surface in clock-wise direction (the diamagnetic 
current) . Depending on which of these currents prevail, 
the magnetisation (or, equivalently, the magnetic mo- 
ment M\ = (l/2c) f\j s r]dv) can change sign, as function 
of h\ (see [8] for details). Recall, that in the vortex-free 
state (m = 0) there exists only diamagnetic current, i.e. 



M\ < 0, see Fig. 1(c). In presence of the vortex (to = 1), 
but in absensce of the external field (h\ = 0), the screen- 
ing current and the magnetic moment correspond to the 
paramagnetic state. This state is metastable, because 
the vortex-free state posesses smaller free-energy, than 
the state to = 1 [2-8]. 

The curve Pq in Fig. 3(a) corresponds to the minimal 
radius R\, when the paramagnetic vortex state (m = 1) 
can still exist inside the homogenious cylinder in absence 
of the field (h\ = 0). [The metastable vortex is held 
inside by the pinning to the cylinder boundary.] In those 
points (i?A,x), which lay below the curve P , to held 
the vortex inside the cylinder it is necessary to impose 
a finite external field, h\ > 0. (This corresponds to the 
field stimulated and re-entrant superconductivity [2-8].) 
Notice, that if the cylinder radius R and the parameter 
x are fixed, to cross the paramagnetic pinning boundary 
(Po) h is sufficient to vary only the sample temperature, 
because R x = R/\(T). 

The presence of a smooth tail in the function ip max {h\) 
[Figs. 3(6) and 4(a)] allows (as in the case to = 0) to use 
the linear approximation (^i < 1) for finding the upper 
boundary of the superconducting state, h c . In the region 
of I-order jumps [s\ in Fig. 3(a)], where the function ip(p) 
is not small, the linear approach fails and more rigorous 
analysis, based on full system of non-linear equations (1)- 
(5), is neccessary. [The boundaries Si_n(x) and C„ s (x) 
themselves can not be found from the linear equation, 
because the latter does not depend on x [11]. The com- 
parison of the results of the rigorous and linear analysis 
will be reported elsewhere.] 

Similarly, one can consider the higher giant-vortex 
states (to > 1, see Fig. 5 for to = 2). Here also ex- 
ist the boundaries of I- and Il-order phase transitions, 
the jumps on the magnetisation curves, the paramag- 
netic and diamagnetic currents, and other peculiarities, 
which are analogous to those presented in Figs. 1-4. 

IV. CONCLUSION AND DISCUSSION 

Basing on self-consistent solution of non-linear system 
of GL-equations, the boundary, Si_n, is found, which 
separates the regions, where the superconducting state of 
the cylinder is destroyed by the external magnetic field 
either by I-order jump (the region s\), or gradually, by Il- 
order phase transition (the regioin sn). This boundary is 
a complicated function of the parameters (to, R\, x) [see 
Figs. 1(a), 3(a), 5]. 

Note, that in the case of an infinite (open) supercon- 
ductor the phase boundary between I- and Il-order tran- 
sitions lays at the value x = 1/ y/2 [9] . (At this value the 
surface energy at the interface of superconducting and 
normal metals vanishes, and the magnetisation M(H) 
acquires a smooth tail [9].) However, the case of infi- 
nite superconductor is degenerated, in the sense that the 
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total number of vortices in the open system can not be 
defined. Due to this degeneration there are many solu- 
tions of the system (1), (2) with different to, and it is 
possible to consider the superconducting state as a linear 
combination of states with different vorticities to [9] . In 
the bounded system this degeneracy is removed, and it 
is necessary to consider the states of fixed vorticity to 
(the quantum number to is now a topological invariant). 
[It is easy to prove that self-consistent solution of Eqs. 
(l)-(5) with ip ^ is unique for every value of h\.} The 
mentiond difference of the si_n-boundary from the value 
x = l/v2 is due to the difference in geometries and to 
the account of the space-quantization effects, present in 
the bounded system. [To trace the limiting transition 
from the bounded to open geometry, it might be necces- 
sary to consider the case of flattend elliptical cylinder, 
which models the geometry of an infinite slab, adopted 
in [9].] _ 

Mention in conclusion, that the main attention in the 
present work was paid to the mathematical side of the 
problem: to describe the si_n-boundary basing on for- 
mal solutions of GL-equations. The important physical 
question of compairing the Gibbs free energies of various 
mathematically possible states, and finding the most sta- 
ble ground state (which the system would occupy in equi- 
librium), was put aside. (Some illustrations of the Gibbs 
free energy behavior, found from Eq. (7), are given, for 
instance, in [6-8].) In justification, it may be reminded, 
that the physical system may occupy not only the groung 
state, but also the excited metastable states of higher 
energy. [In particular, the controversial paramagnetic 
Meissner effect may be attributed to the metastable vor- 
tex states in the mesoscopic system, see [8] for details). 
The formal solutions of GL-equations describe all possi- 
ble states, including stable and metastable one, thus the 
full analysis, based on these equations, may be pertinent 
to the experiment. However, the case of infinitely long 
cylinder, considered above, approximates rather poorly 
the geometry, used in real experiments. In this respect, 
the superconducting disks, considered in [4,5], are more 
adequate, though it would be more difficult to obtain 
rigorous solutions for the bounded 3-dimensional sam- 
ple. Thus, further analysis of the questions, touched on 
in the present paper (as well as possible connection with 
experiment), is necessary. 
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Figures captions 

Fig. 1. (a) - The boundary (si_n) between regions (si 
and sn), where I- or Il-order phase transitions to normal 
state (V> = 0) in magnetic field occure. (b) - The behavior 
^max{h\) in the points 1-6 (to = 0, R\ = 4) in Fig. 1(a). 
In the region si the order parameter vanishes by I-order 
jump. In the region sn the order parameter ipmax(h\) 
has a "tail", and vanishes smoothly, by Il-order phase 
transition, (c) - Analogous behavior for magnetisation, 
M\(h\). The peep-holes 1-9 in (a) are pearced in the 
points: l-x = 0.2, 2-X = 0.4, 3-x = 0.7, 4 - x = 1, 
5 - x = 1.0.5, 6 - X = 1.2 (R x = 4); 7 - R x = 3, 8 - 
RX = 2, 9- RX= 1.5 (x= 0.7). 

Fig. 2. The dependencies: (a) - ipmax(h\) and (b) - 
M\(h\) for m = 0, x = 0.7. The numeration of curves 
correspond to points 3,7-9 in Fig. 1(a). 

Fig. 3. Analogous to Fig. 1, but for m = 1. The 
dashed curve C ns in Fig. 3(a) separates the normal (n) 
and superconducting (s) regions. The curve P$ marks 
the points (R\,x), where the metastable vortex state 
(to = 1) may still exist in absence of the field (h\ = 0) 
due to the pinning to the boundary. Below the curve 
Pq the vortex state may exist only in presence of finite 
external field (h\ > 0, see the curves 1,2 in Figs. 3(o, c)). 
(This is an example of the field stimulation effect, or re- 
entrant superconductivity.) The peep-holes 1-8 in (a) are 
pearced in the points: 1 - x — 0.35, 2 — x = 0.4, 3 - 



4 



x = 0.7, 4 - x=l, 5- x= 1.07, 6 - h = 1.2 (i? A = 4); 
7- i? A = 3, 8- R x = 2.4 (x = 0.7). 

Fig. 4. Analogous to Fig 2, but for m = 1. The 
presence of paramagnetic (M A > 0) and diamagnetic 
(M\ < 0) parts of magnetisation is evident in Fig. 4(6). 

Fig. 5. Analogous to Figs. 1(a) and 3(a), but for 
m = 2. The vertical asymptote x = 0.94 is the same 
for m = 0,1,2. This is natural, because for large radii 
(Rx 3> 1) the influence of the vortex field is negligible. 
The bottom of the curve Si_n lays at Rx = 2.78 (with 
Rx = 2.45 for m = 1, and R x = 1.69 for m = 0). The 
dashed curve C„ s is well approximated by the dependence 
C ns w 1.81/x (the dotted line). 
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